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Abstract 



In this chapter we summarize recent developments in the study of kinetically con- 
strained models (KCMs) as models for glass formers. After recalling the definition of 
the KCMs which we cover we study the possible occurrence of ergodicity breaking 
transitions and discuss in some detail how, before any such transition occurs, relax- 
ation timescales depend on the relevant control parameter (density or temperature). 
Then we turn to the main issue: the prediction of KCMs for dynamical heterogeneities. 
We focus in particular on multipoint correlation functions and susceptibilities, and de- 
coupling in the transport coefficients. Finally wc discuss the recent view of KCMs as 
being at first order coexistence between an active and an inactive space-time phase. 
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1.1 Motivation 

Kinetically constrained models (KCMs) are simple lattice models of glasses. They fur- 
nish a perspective on the glass transition problem which has its origin in the work of 
Glarum ( |Glarum, 1960D , Anderson( Anderson, 1979 1 and coworkers ( |Palmer et al., 1984| , 



and Andersen and coworkers ( |Fredrickson and Andersen, 1984[ ). This perspective as- 
sumes that most of the interesting properties of glass forming systems are dynam- 
ical in origin, while thermodynamics plays a very limited role. KCMs tend to have 
simple and uninteresting thermodynamics, typically that of a non-interacting lattice 
gas. In contrast, they display rich dynamical behaviour as a consequence of kinetic 
constraints. This combination of simple thermodynamics and locally constrained dy- 
namics is often assumed to be the result of coarse-graining of a dense molecular 
system flGarrahan and Chandler, 2003 ): dense fluids are structureless at distances be- 



yond the molecular length, but interatomic forces at high densities are highly con- 
straining, giving rise to local restrictions in the dynamics. As such KCMs are meant 
as models of glass forming systems at high densities or low temperatures, and aim 
to capture their dynamical behaviour for motion beyond the inter-molecular distance 
and for long times. 

KCMs generally use local constraints. Nevertheless, as we discuss below, these give 
rise to collective dynamics due to a form of dynamical frustration: at low tempera- 
tures/high densities there is a conflict between the scarcity of excitations/vacancies 
and the need for them to "facilitate" local motion, leading to hierarchical and cooper- 
ative relaxation. Note that this frustration, in contrast with e.g. the random first-order 



approach ( |Lubchenko and Wolynes, 2007[ ), does not arise from quenched disorder, an 
ingredient which at any rate is not trivial to justify in models of real liquids. A fur- 
ther distinction from the latter approach is that KCMs offer a "non-topographic" 
( |Berthier and Garrahan, 2003[ ) view of the glass transition problem. By this we mean 
that it is not a change in the topographic structure of the potential energy landscape 
(such as a transition between a saddle-dominated and a minima-dominated regime) 
that determines glassy relaxation, but a change in the degree of connectivity of the 
configuration space. Indeed, due to the presence of the constraints, the effective con- 
nectivity of this space depends on temperature/density, as we discuss in more de- 
tail in section 11.31 This picture provided by the KCM approach is appealing as an 
effective description of the physics for example in the case of hard spheres, where 
all allowed configurations have the same energy. For a more detailed discussion on 
the non-topographic KCM approach versus the energy landscape scenario (and addi- 



Motivation vii 



tional references on the latter) we refer the reader to flBerthier and Garrahan, 2003| 
|Whitelam and Garrahan, 2004[ ). Further information on disorder or landscape based 
approaches, along with other explanations of glassy behavior that invoke a thermody- 
namic transition, can be found in Chapter 1 of this book. 

KCMs are simple enough to allow for detailed analysis. The current interest in 
KCMs originates in the fact that they exhibit explicit mechanisms for super- Arrhenius 
slowdown and stretched relaxation ( |Sollich and Evans, 1999[ ) as a consequence of lo- 
cal, disorder free interactions, and without the emergence of finite temperature sin- 
gularities (Toninelli et al., 20046). At the same time they provide a natural explana- 
tion (Garrahan and Chandler, 2002) for the phenomenon of dynamical heterogeneity 

(for reviews see ( |Ediger, 2000| |Glotzer, 2000| |Andersen, 2005|». KCMs are to the con- 

strained dynamics/facilitation perspective ( Palmer et al., 19 84; Fredrick son and Andersen, 1984[ ) 
what the random energy model and the p-spin spin glass are to the random-first or- 
der transition approach ( |Lubchenko and Wolynes, 2007[ ) . Insights from the analysis of 
KCMs allow one to construct a comprehensive theoretical picture of the glass transition 
problem (see (Chandler and Garrahan, 2010) for a review) which is quite distinct from 
other competing theories flGotze and Sjogren, 199 2 Debc nedetti and Stillinger, 2001| 
|Lubchenko and Wolynes, 2 007 Kivclso n and Tarjus, 2008: Cavagna, 2009). Beyond their 
importance as models for describing glassy phenomenology, KCMs are of interest for 
the mathematical community. Indeed, even though they belong to the class of interact- 
ing particle systems with Glauber and Kawasaki dynamics, the rate at which elemen- 
tary moves (birth/death or jump of particles) occur may degenerate to zero due to the 
presence of the kinetic constraints. This prevents the use of the standard probabilis- 
tic tools developed for such systems. Furthermore it gives rise to peculiar phenomena 
including the presence of several invariant measures, ergodicity breaking transitions 



(Toninelli and Biroli, 2008a), unusually long mixing times (Cancrini et al, 2008) and 



aging phenomena (Faggion ato et al, [ ). 

The most recent comprehensive review of KCMs is that of Ritort and Sollich 
dRitort and Sollich, 20011 ) which covers the literature until roughly the end of 2001. 
Earlier surveys can be found in Refs. flJackle, 1986| |Fredrickson, 1988] |Palmer, 1989[ ). 
The aim of this chapter is to briefly describe the developments in the last decade or 
so. It is organised as follows. In Section H~2l we introduce the KCMs that we will cover, 
grouping them according to whether they are non-conservative (Glauber dynamics) or 
conservative (Kawasaki dynamics) . Section 11.31 is devoted to a discussion of probably 
the most basic mathematical manifestation of glassiness in KCMs, i.e. whether they 
possess ergodicity breaking transitions. Beyond such a transition, the relaxation time 
to equilibrium is effectively infinite. Physically, of course, this may not be distinguish- 
able from a finite but very long time, and so we consider in Section [1.41 how relaxation 
timescales in KCMs depend on the relevant control parameters (temperature/density). 
These times do indeed generically increase extremely fast, in particular in the so-called 
cooperative models. In Section 11.51 finally, we turn to the predictions of KCMs for 
dynamical heterogeneity in glasses. Here the models provide a appealingly intuitive 
picture that is directly based on considering the dynamics in real (rather than e.g. 
Fourier mode) space. We cover definition and predictions for multi-point correlation 
functions and susceptibilities, decoupling in the temperature or density-dependence of 
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transport coefficients, and finally the view of KCMs as being at a first-order coexis- 
tence between two different dynamical (space-time) phases. We conclude in Section fTTBl 
with a summary and a critical discussion of the advantages and drawbacks of KCMs. 
We also provide a brief comparison with alternative approaches, and finally an outlook 
towards future work. 



1.2 The models 

The majority of KCMs are defined as stochastic lattice models with binary degrees of 
freedom, and it is on such models that we focus in this chapter. We do not attempt to 
give here an exhaustive survey of the wider range of models studied previously, but refer 



to this for the review listed in the introduction (Ritort and Sollich, 2003 Jackie, 1986 
|Fredrickson, 1988] |Palmer, 1989[ ). For our present purposes, then, a KCM has on each 
lattice site i an occupation variable rii € {0, 1}, and the collection of the rii defines the 
overall configuration n. Conservative KCMs are lattice gases where the rii indicate the 
presence (rii = 1) or absence (rii = 0) of particles. The dynamics follows a continuous- 
time Markov process that consists of a sequence of particle jumps, subject to kinetic 
constraints as explained below. The total particle number rii is conserved. 

We will also discuss non-conservative KCMs. These are motivated by a conceptual 
coarse-graining to a length scale several times the particle diameter of the underlying 
physical glass that is being modelled. Each lattice site i then represents a small region 
of material containing at least a few particles, and one sets n,- = 1 (respectively n, = 0) 
if the density in this region is above (respectively below) a certain threshold where 
local re-arrangements - of the particles inside the element - become possible. This 
representation no longer contains the precise values of all local densities, and conse- 
quently in non-conservative KCMs moves that change J^. n.; are allowed. The reader 
should be aware that in much of the literature a reverse convention for the two states 
of rii is used, with n.; = 1 standing for high mobility, i.e. low density, and vice versa 
for ni = 0. (These states are then often further identified as up- and down-spins, but 
we will avoid this terminology.) We opt for the version discussed above as it makes for 
a unified discussion of conservative and non-conservative KCMs: = 0, an empty or 
low-density site, is mobile in both contexts and facilitates the dynamics on neighbour- 
ing sites. For brevity we use the terms "empty" and "low-density" interchangeably 
below, and similarly for "occupied" and "high-density" . When we talk about the den- 
sity of the system, we correspondingly mean the fraction of sites with m = 1, both in 
conservative and non-conservative KCMs. 

The key property of all KCMs is that in order to perform a move the configura- 
tion must satisfy a local constraint which usually corresponds to requiring a minimal 
number of empty sites in some appropriate neighbourhood. This represents the phys- 
ical intuition that, when particles re-arrange in a glass, motion in any given region 
requires the presence of mobile regions around it ( |Fredrickson and Andersen, 1984| 
|Kob and Andersen, 1993] ). When the density of these mobile or facilitating sites de- 
creases, the dynamics slows down. Another property shared by all models is that the 
rates satisfy detailed balance with respect to (w.r.t.) a Boltzmann distribution that 
factorizes over sites (mathematically, a Bernoulli product measure), so that there are 
no static interactions, beyond the effective hard core repulsion implemented by the 
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restriction rii € {0, 1}. As explained in the introduction, this idealization means that 
KCMs can also be viewed as an attempt to find out how much of glass phenomenology 
can be explained on purely dynamical grounds, without recourse to e.g. static phase 
transitions. 

KCMs as described above are plainly quite simplistic compared to more realistic 
interacting atomic or molecular systems. As we shall explain, however, for appropriate 
choices of the constraints they display a behavior which is in agreement with the 
broad phenomenology of glass forming liquids including super- Arrhcnius slowing down 
of the dynamics, stretched exponential relaxation, dynamical heterogeneities, aging 
phenomena and ergodicity breaking transitions. 



1.2.1 Facilitated spin models: FA, East and Spiral models 

We next describe the non-conservative KCMs that we will consider. Because mobile 
(rij = 0) sites facilitate motion on neighbouring sites, and because binary degrees of 
freedom with non-conservative dynamics can be interpreted as spins, such models are 
also called "facilitated spin models" . In these models, the rate for changing the state 
of site i is / 2 ;(n)[(l — p)ni + p(l — rii)}, where fi depends on the configuration n in 
a finite neighbourhood of i, but not on rii itself. It is then immediate to verify that 
detailed balance holds w.r.t. a Boltzmann distribution with energy function — J^. m 
and inverse temperature (3 linked to the density by p = 1/(1 + e^* 3 ): in equilibrium we 
have rii = 1 with probability p and n% = with probability 1 — p, independently at 
each site. In general fi is chosen to be nonzero (allowed move) when an appropriate 
neighbourhood of i contains a minimal number of empty sites, otherwise it is zero 
(forbidden move). As the temperature decreases, the density p of occupied sites 
increases and the density q = 1 — p = 1/(1 + e* 3 ) of empty sites decreases: the dynamics 
must then slow down. 

Non-conservative KCMs can be divided into two classes: non-cooperative and coop- 
erative models. For the former it is possible to construct an allowed path - a sequence 
of configurations linked by transitions with nonzero rates - which completely empties 
any configuration provided that it contains somewhere an appropriate finite cluster of 
empty sites, the mobile defect. For cooperative models, no such finite mobile defects 
exist. As we will detail in Section fl.4| non-cooperative and cooperative models display 
respectively an Arrhcnius and super- Arrhenius slowing down. 

We now give the definitions of some specific models which cover these two differ- 
ent classes of behaviour. The first class of non-conservative KCMs was proposed by 
Fredrickson and Andersen (FA) ( |Fredrickson and Andersen, 1984 1, hence the name 



Fredrickson- Andersen models. The simplest of these models, which we write as FA-1 
("one-spin facilitated FA"), allows a change of state at site i only if at least one of 
the nearest neighbours is empty: fi(n) = 1 if X^~;(l — n i) > 0, fi{n) = otherwise, 
where the sum runs over the nearest neighbours j of site i. It is easy to check that 
the presence of a single empty site allows one to empty the whole lattice: the model 
is non-cooperative. For theoretical calculations it is often easier to define fi somewhat 
differently, as fi{n) = X^~i(l — n j)- This leaves the kinetic constraints enforced by 
vanishing rates as they are and just changes the rates for allowed moves when more 
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than one empty neighbour is present; as such, it makes no qualitative difference to the 
behaviour of the model. 

One can similarly define m-constrained FA models, FA- to, by setting fi(n) = 
1 if at least to neighbouring sites of i are empty, and fi(n) = otherwise. These 
models are normally considered on hyper-cubic lattices of dimension d, with 2 < to < 
d flFrcdrickso n and Andersen, 1984[ ). As can be checked directly, it is not possible in 
any of these models to devise a finite cluster of empty sites which always lets one empty 
the entire lattice. Consider for example the case d = 2, to = 2 in infinite volume and 
focus on a configuration which contains two adjacent infinite rows of occupied sites. 
Even if the rest of the lattice is completely empty, none of the sites in these two 
rows can ever change state because each has at most one empty neighbour. Thus it is 
not possible to devise a mobile defect which can unblock (empty) every configuration 
and the model is cooperative. The same conclusion applies on a finite lattice, e.g. with 
periodic boundary conditions. The restriction on to comes from the fact that, if to > d, 
it is possible to construct finite structures of occupied sites that are blocked, i.e. can 
never change state. They are therefore not suitable to describe the slow dynamics close 
to a glass or jamming transition because a finite fraction of the system is jammed at 
any density. 



Another example of a cooperative model is the one-dimensional East model (Eisinger and Jackie, 1991 1 
In this case the constraint requires for an allowed move that the left nearest neighbour 
be empty, that is /; = 1 — rii— \. Note that on a finite chain the presence of a single 
empty site at the left boundary lets one empty the entire chain; but this is not a mobile 
defect because it will work only when it is found in a specific position. In a chain that 
is full except for a finite cluster of empty sites in a generic location, the part of the 
chain to the left of the cluster cannot be emptied, and so the model is cooperative. 
One can show that due to the directed nature of the constraint the relaxation involves 
the cooperative rearrangements of increasingly large regions as q becomes small. This 
leads to super- Arrhenius behaviour of relaxation time scales ( [Evans and Sollich, 1999| 
Sollich and Evans, 2003[ ) . (The East model shows up a slight shortcoming of our def- 



inition of co-operativity: with periodic boundary conditions, we would formally label 
this model non-cooperative because all sites can be emptied starting from a single 
empty site. Physically, the model is clearly co-operative independently of the bound- 
ary conditions.) We discuss below also higher-dimensional generalizations of the East 
model, e.g. in d = 3 the North-East-Front (NEF) model, where a local change of state 
is allowed if at least one of the neighbouring sites to the North or East or front is 
empty. 

The models described above have been closely studied since they were proposed in 
the 80s fFredrick son and Andersen, 1984[ ) and 90s ( Eisinger and Jackie, 199lj ), with 
some work continuing into the present decade ( |Ritort and Sollich, 2003 1, particu- 
larly regarding the non-equilibrium behaviour (see Section 11.4.3)) and detailed com- 
parisons with glass phenomenology (Garrahan and Chandler, 2003 ] | Jung et al., 2004 



|Berthier et al, 2005a|[Toninelli et al, 20056[|Chandler et al, 2006||Elmatad et al, 2009D 
In parallel, progress has come from the development of new KCMs dGarrahan and Chandler, 2003[ ). 
We discuss below one of these, which reproduces in particular the mixed character of 
the glass transition, where diverging lengthscales characteristic of second order tran- 



The models xi 



sitions combine with a first order jump in the fraction of frozen degrees of freedom. 
This two-dimensional Spiral model ( |Toninelli et at, 2007| |Toninelli and Biroli, 2008~a| 
Toninelli and Biroli, 200861 is defined as follows. Consider a square lattice and, for 



each site i, divide its first and second neighbours into the North-East (NE), South- 
West (SW), North- West (NW) and South-East (SE) pairs, as illustrated in Fig. Ol 
Then the constraint for a move to occur at i is that "both its NE and/or both its 
SW neighbours should be empty" and "both its SE and/ or both its NW neighbours 
should be empty" (see Fig. in other words, fi(n) = 1 if this condition is sat- 

isfied and = otherwise. Stated in a more geometrical form, the constraint requires 
that all four sites in one of four contiguous sets of neighbours (NE and SE, or SE 
and SW, or SW and NW, or NW and NE) have to be empty. The motivation for 
the somewhat involved form of this constraint will become clear when we analyse the 
dynamical transition in the model below. We note that before the Spiral model a 
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Fig. 1.1 Site i and its NE, NW, SE and SW neighbours. The constraint is (is not) verified 
at the central site in case (a) (in case (b)). 



somewhat more involved choice of the constraints (the so called Knights model) was 
proposed ( |Toninelli et at, 2006[ ?) which leads to the same mixed first /second order 
transition. Other choices in the same spirit have also been investigated since, see for 



example (Jeng and Schwarz, 



1.2.2 Kinetically constrained lattice gases: KA and TLG models 

Conservative KCMs are defined in a similar spirit to their non-conservative analogues, 
but with dynamics that conserves total particle number. A particle at i attempts at a 
fixed rate (= 1 without loss of generality) to jump to any empty nearest neighbour site 
j. The rate for this move is thus rii(l — rij) (n) . The factor fij(n) again implements a 
kinetic constraint and is chosen not to depend on the configuration at sites i and j. The 
resulting dynamics preserves the number of particles, and so detailed balance on finite 
lattices is satisfied w.r.t. the Boltzmann distribution which is uniform on configurations 
with the appropriate fixed particle number. By linearly combining distributions with 
different particle numbers one can of course also create grand-canonical Boltzmann 
distributions where each site independently contains a particle with some probability 
p. Detailed balance is also satisfied w.r.t. these distributions, and they extend naturally 
to the limit of an infinite lattice. 
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Again, we can classify conservative KCMs into non-cooperative and cooperative 
models. For the former it is possible to construct a finite group of empty sites, the 
macrovacancy, such that for any configuration the macrovacancy can be moved all over 
the lattice and any jump of a particle to a neighbouring empty site can be performed 
when the particle is adjacent to the macrovacancy. Clearly the baseline lattice gas 
in which there are no kinetic constraints, which in d = 1 is the symmetric simple 
exclusion process (SSEP), is non-cooperative and the minimal macrovacancies are just 
isolated empty sites. 

Chronologically the first conservative KCM is the Kob-Andersen model (KA) 
( |Kob and Andersen, 1993| , which has cooperative dynamics. Here a particle can jump 
to a neighbouring site only if both in the initial and final position at least to of its 
nearest neighbouring sites are empty. The original KA model had m — 3, with par- 
ticles on a cubic lattice, but one can similarly define KA-m models on hypcr-cubic 
lattices of dimension d and for different values of the parameter to, with 2 < m < d. 
The restrictions on to arise from the fact that to = 1 corresponds to an unconstrained 
lattice gas, while a model with to > d has a finite fraction of frozen particles at any 
density (as for FA models in the same parameter regime). Similarly to the case of 
FA-TO models with m > 2, one can directly check that all KA-m models are coop- 
erative. On other lattices, KA models can be non-cooperative. Consider for example 
the KA model with m = 2 on a triangular lattice. Here one can check that a "dimer" 
of two neighbouring empty sites forms the required macrovacancy, which can move 
across the lattice in a tumbling motion even when all other sites are occupied. Two 
other KCMs on triangular lattices w ere introduced in (|Jackle and Kronig, 1994[ ) and 
have more recently been analysed in (Hedges and Garrahan, 2007). For the (l)-TLG a 
particle can move from site i to a neighbouring site j if at least one of the two mutual 
neighbours of i and j is empty. For the (2)-TLG a particle can move from i to j if 
both the two mutual neighbour sites are empty. It is easy to verify that the (l)-TLG 
is non-cooperative, with - as for KA with m = 2 on the same lattice - a dimer of 
vacancies being a macrovacancy. The (2)-TLG, on the other hand, is cooperative: a 
chain of particles occupying an entire row of the lattice, for example, can never be 
destroyed, thus it is not possible to construct a macrovacancy. 



1.3 Ergodicity breaking transitions 

As explained above, the equilibrium distribution in KCMs is trivial because it factorizes 
over sites. In particular, then, no equilibrium phase transition can occur. On the other 
hand, the presence of constraints might induce transitions of purely dynamical type: 
detailed balance alone does not guarantee that the distribution over configurations will 
converge for large times to the Boltzmann equilibrium distribution. To see this, return 
to the example of the FA-2 model on a square lattice (m = d = 2), either finite with 
periodic boundary conditions or infinite, and a starting configuration n(0) which has 
two adjacent rows of occupied sites. Then since these sites are forever blocked, for any 
site i within the two rows, lim t _ >00 (ni) t = 1 ^ (n,)eq = p, where (-)t, is the average over 
the stochastic dynamics which starts from n(0) and (-) eq is the equilibrium average. 

Motivated by this, wc can ask whether convergence to the equilibrium distribution 
is recovered at least in the thermodynamic limit: if we sample an initial configuration 
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n(0) from the equilibrium distribution, then is the large time limit of the average 
of any function g(n) equal (with probability one) to its equilibrium average, i.e. is 
\im t ^oo(g)t = (.g)cq? If this is the case we will say that the system is ergodic. The 
models we consider are all ergodic at sufficiently high density q = 1 — p of facilitating 
(empty) sites. If they do become non-ergodic as q is reduced, we call q c the critical 
density of empty sites at which this transition occurs. As will be discussed below, the 
models defined in the previous sections are ergodic at any positive q = 1 — p, so that 
q c = 0, with the exception of the Spiral model and FA models on Bethc lattices which 
display ergodicity breaking transitions at < q c < 1. 

Owing to the factorized form of the equilibrium distribution, ergodicity corre- 
sponds to the fact that in the thermodynamic limit the configuration space is covered 
by a single irreducible set, i.e. a set of configurations which are connected to each 
other by allowed paths; for a proof see ( |Cancrini et al., 2008[ Proposition 2.4) and 
( |Cancrini et al, 2010a| Section 2.3). This in turn is equivalent to the requirement that 
for any site i (in non-conservative KCMs) or pair of nearest neighbour sites i,j (in 
conservative KCMs) there is an allowed path which transforms the configuration into 
one where the constraint is satisfied at i In other words, the probability that any 

site i belongs to a cluster of forever blocked sites must vanish in the thermodynamic 
limit. 

For non-cooperative KCMs one sees easily that ergodicity holds at any density 
p < 1, i.e. at any q > 0. In the non-conservative case, the probability of finding at 
least one mobile defect in an equilibrium configuration goes to one in the infinite 
volume limit, and then starting from this defect one can empty all sites. Thus, any 
configuration is connected to the "all empty" configuration with probability one. The 
situation is similar for non-cooperative but conservative KCMs. 

The case of cooperative KCMs is more delicate and here an ergodicity breaking 
transition can occur. The non-conservative case is again simpler, so we begin with 
this. In order to determine the probability that a site belongs to a blocked cluster, 
consider the following deterministic procedure: iteratively empty all sites for which 
the constraint is verified until we reach either the completely empty configuration, or 
one in which there is a "backbone" of mutually blocked occupied sites. Then it is easy 
to verify that the sites that are blocked forever under the stochastic dynamics of the 
actual KCM are precisely those that are empty at the end of this deterministic proce- 
dure. In other words, the problem of the existence of an ergodicity breaking transition 
for cooperative non-conservative KCMs can be reformulated as a percolation transition 
for the final configuration of the above deterministic dynamics. For FA-m models this 
deterministic dynamics coincides with the well known algorithm of bootstrap perco- 
lation and the results in (Aizenman and Lebowitz, 1988 Schonmann, 1992) establish 
that q c = on hypcrcubic lattices in any dimension d and for any facilitating pa- 
rameter 1 < m < d (while trivially q c = 1 for m > d since finite blocked structures 
can occur). This work disproved a long-standing conjecture from the original FA pa- 
per (Fredrickson and Andersen, 19841, namely that an ergodicity breaking transition 
would occur at some q c > 0. Such a transition does take place, however, when one 
considers FA models on a Bethe lattice as will be explained in Section ll.3.11 An exam- 
ple of a finite-dimensional model displaying such a transition is provided by the Spiral 
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model (see Section ["l.3.2p . 

For cooperative conservative KCMs, the proof of ergodicity is more involved. For 
example, for the KA-2 model on a square lattice (d = 2), we have shown that 
the irreducible set of configurations that has unit probability in the thermodynamic 
limit is the one containing all configurations which can be connected by an allowed 
path to a configuration which has a frame of empty sites on the last shell before 
the boundary ( |Toninelli et al, 2005 a| ). Establishing that this set has unit probabil- 
ity in the thermodynamic limit is more complicated than in the corresponding non- 
conservative KCMs, i.e. FA models. This is because there is no deterministic bootstrap- 
like procedure which allows one to establish whether or not a configuration does 
or does not belong to the irreducible set. Nevertheless, a formal proof can be con- 
structed ( |Toninelli et al, 200 5a; Can crini et al, 2010aj ), and demonstrates that the er- 
godicity breaking transition originally conjectured by Kob and Andersen flKob and Andersen, 1 993) 
does not exist. Analogous arguments can be constructed for the other choices of m 



and d (Toninelli et al, 2005a| to rule out the occurrence of a transition. 



1.3.1 FA models on Bethe lattices 

As a simple example of a model that does have an ergodicity breaking transition we 
consider next the FA-m model on a Bethe lattice, i.e. a random regular graph of 
connectivity k + 1. Exploiting the local tree-like structure of such a graph, it is easy 
to write a self-consistent equation for the probability P that a given site is occupied 
and blocked, conditional on the fact that its ancestor is occupied and blocked: 

rn-l , k 



p = (i-?)]T p^i-p)*. (i.i) 

i=0 

The factor 1 — q is the probability for the given site to be occupied; the sum gives the 
probability that at most m — 1 of the descendants of the site are empty so that it is in 
fact blocked. The leading term on the r.h.s. for small P is 0(P k ~ m+1 ). For m = k, the 
largest value that does not produce finite blocked structures, this is linear and so one 
gets a continuous transition at (1 — q c )m = 1. For all smaller m, the small P-increase 
is with a higher power of P and so a discontinuous transition results, with P = for 
q > q c and P = P c + 0(y/q c — q) for q < q c , where P c > 0. The mechanism behind 
this combination of a discontinuous onset and a critical singularity has been anal- 



ysed in detail (S N Dorogovtsev and Mendes, 2006 J M Schwarz and Chayes, 2006) 



In particular, the singular square root behavior is due to the extreme fragility of the 
infinite spanning frozen or "jammed" cluster at the transition, and to the existence of 
a related length scale which diverges as a power law. 

1.3.2 Spiral model 

In this section we will explain the mechanism behind the ergodicity breaking transition 



of the Spiral model flToninelli et al, 2007| ?; |Toninelli and Biroli, 2008 § ). Consider the 
directed lattice that is obtained from the square lattice by putting two arrows from each 
site towards the neighbours in the NE pair, i.e. pointing North and North-East. The 
resulting lattice is a tilted and squeezed version of a two-dimensional oriented square 



Ergodicity breaking transitions xv 



lattice (see Figure [L2|l . Therefore, if the density is larger than the critical density 
of oriented site percolation (directed percolation, DP), pj? p , there exists a cluster of 
occupied sites which spans the lattice following the direction of the arrows. Consider 
now a site in the interior of this directed cluster (see Figure [O]): by definition there 
is at least one occupied site in both its NE and SW neighbouring pairs, therefore the 
site is blocked with respect to the constraints of the Spiral model. Thus for p > pj? p 
the system contains one or more blocked clusters and is therefore non-ergodic. This 
suggests that the ergodicity breaking transition occurs at q c = 1 — pj? p , but the 
argument so far does not exclude a transition earlier, i.e. at a larger density q of empty 
sites. Indeed, since blocking can occur along either the NE-SW or the NW-SE direction 
(or both) , the presence of a blocked cluster does not imply that there is a directed path 
through the lattice as in the DP argument. To establish that, nevertheless, blocked 
clusters do not occur for q > 1 — pj? p , one shows that empty regions of linear size 
much larger than the parallel length of DP, ^ , act as "critical defects" : starting from 
any such empty region we can very likely empty the whole lattice. Because £|| is finite, 
so is the size of these defects. In the thermodynamic limit at least one such defect 
will occur in the system, and so for q > 1 — p® p any configuration can be completely 




Fig. 1.2 The tilted and squeezed two dimensional oriented lattice obtained by drawing 
arrows from each site to its neighbours to the North and East. The site marked by the cross 
belongs to an oriented occupied cluster and is therefore blocked. 

An important feature of the ergodicity breaking transition in the Spiral model is 
that it is discontinuous, in the sense that the fraction of blocked sites jumps to a 
nonzero value at q = q c . This is relevant for the connection to real glasses, because 
it implies that two-point correlation and persistence functions must display a plateau 
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in their time-dependence when the transition is approached (q > q c ), as is observed 
experimentally. Earlier models with a transition at a finite q c , e.g. a two-dimensional 
generalization of the East model ( |Reiter et ai, 1992] ) , develop only a fractal cluster of 
blocked particles at the transition, which occupies a vanishing fraction of the system. 
The fact that in the Spiral model blocked clusters are compact at q c follows by a direct 
construction of blocked structures ( |Toninelli et ai, 2007| |Toninelli and Biroli, 2008 a| 
|Toninelli and Biroli, 20086| and is a consequence of the presence of the two transverse 
blocking directions in the constraints, and of the anisotropy of DP. 

1.3.3 Summary: Presence/absence of ergodicity breaking transition 

For the sake of clarity let us summarize the results of this section by listing the KCMs 
we have discussed according to whether or not they display an ergodicity braking 
transition at a non-trivial (different from or 1) critical defect density q c . 

The models that do display such a transition are the FA and KA models on Bethe 
lattices, the Spiral model and the North-East model. The latter was not mentioned 
above because it is not of direct interest for modelling the glass transition, due to the 
continuous character of its ergodicity breaking transition. It is a spin facilitated model 
on a square lattice with the constraint requiring both the North and East neighbour to 
be empty. It is easy to verify along the lines of the argument in the previous subsection 
that the North-East model as defined by this constraint displays a transition at the 
critical density of oriented percolation. 

The models that do not display a transition arc: the FA-m and KA-m models for 
any choice of m and on hypercubic lattices in any finite dimension d; the East model; 
and the (l)-TLG and (2)-TLG models. 



1.4 Bulk dynamics of KCMs 

1.4.1 Glassy time scale divergences and static length scales 

For KCMs to be useful as models of physical glasses they need to satisfy the basic 
requirement of dynamical slowing down as the glassy regime approached. In our case 
this corresponds to increasing density p, or decreasing density of empty sites q. We will 
therefore now give an overview of numerical and analytical results on how relaxation 
time scales r diverge in KCMs when q approaches q c . In concrete terms, we will take 
t as the typical time in the relaxation of density-density correlation and persistence 
functions. In cooperative KCMs, r turns out to be connected to statically defined 
blocking lengths, as we also discuss. 

Starting with FA-1 models, relaxation occurs via the effective diffusion of empty 
sites: an empty site facilitates the emptying of a neighbour site, with rate q = 1— p, and 
the original site can then fill and does so with probability 1/2 before the new empty 
site does. Empty sites can thus diffuse freely, with an effective diffusion constant q/2. 
This suggests that r should grow as an inverse power of q. Because q ~ exp(— /?) 
at low temperatures 1//3, this corresponds to Arrhcnius scaling. Indeed, by an exact 



mapping to a diffusion-limited aggregation model one can derive (Jack et ai, 20066) 
that t ~ l/q z with z = 3 in d = 1 and z = 2 in d > 2. The d = 1 result is 
simple to understand: empty sites are typically a distance 1/q apart, and relaxation 
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requires diffusion across this distance with the effective diffusion constant q/2, so 
T~ {l/ q ?/(q/2) -1/q 3 . 

For FA-m with m > 1, an isolated empty site is unable to move on its own and has 
to wait for a mobile defect, i.e. an appropriate region of empty sites, to move coopera- 
tively into its neighbourhood and so facilitate its motion. As we will detail below, the 
typical number of moves involved in this cooperative process increases as q — > 0, thus 
a super- Arrhenius scaling of the relaxation time has to result. This has been confirmed 
by several numerical investigations which have proposed different forms for the density 
dependence of r ( |S Butler, 1991| |Fredrickson and Brawer, 1986| |Graham et al, 19"97| . 
In order to better understand the cooperative mechanism let us consider for example 
the case m = d = 2. An I x t square of empty sites can be expanded by one lattice 
spacing in both directions, provided that at least one empty site is present on two 
adjacent sides of the square. This is very likely to be the case when t ^> l/q. On the 
other hand, the probability of being able to empty all sites in a region of size l/q with- 
out using external empty sites can be shown to be proportional to exp(— c/q) with c 
a constant of order one ( Aizenman and Lebowitz, 1988; Holroyd, 2003). This result is 
obtained by iterating the process outlined above, multiplying the relevant probabilities 
at each step: one can iteratively remove all the particles in any region starting from its 
interior provided there is a central 2x2 square of empty sites and one additional empty 
site somewhere on each of the four sides of each subsequent shell. In ( jReiter, 1991[ ) 
it was conjectured that relaxation occurs via the diffusion of these critical defects 
leading to a relaxation time diverging as td/pd, with pu — exp(— c/q) the defect 
density and tjj their diffusion time. The latter, as detailed in (Tonincll i et al., 20054 
Section 6), should be a sub- leading correction, giving r ~ 1/pd to leading order. In 
principle one might think that this is not the optimal relaxation mechanism and that 
other relaxation processes could be much faster, avoiding the super- Arrhenius scaling. 
However, the results in (Aizenman and Lebowitz, 1988: Holroyd, 2003) also give the 
typical size of the incipient blocked cluster to which any fixed site belongs and that 
has to be eroded via successive moves from its boundary before the site in qustion 
can be unblocked. This size, L c , diverges as exp(c/2g) thus providing (due to finite 
speed of propagation) a lower bound on time scales of the form t > exp(c/2q) which 
confirms the super- Arrhenius scaling. Note that l? c scales as the inverse of pn, which 
(as the probability that a region of site l/q x l/q is a critical defect) is essentially 
the number density of critical defects. Indeed, clusters of linear size larger than L c 
are typically unblocked because they contain at least one critical defect from which 
relaxation can occur. In (Schonmann, 1992) the typical size of such incipient blocked 
clusters was derived for generic m and d leading to r growing at least as fast as 
exp°( m_1 ) (c/q ,1 /' d-m+1 - ) ). Here exp os is the exponential iterated s times, so that the 
divergence of the time scale is extremely rapid for any m > 3. 

For KA models the basic relaxation mechanism is very similar ( |Toninelli et al., 2005 a[ ) , 
the main difference being that the critical defects arc now regions in which empty sites 
and particles are arranged in such a way that one can find an allowed path - a se- 
quence of allowed moves - to perform any nearest neighbour exchange. As shown in 
dToninelli et al., 2005 a| |Toninelli et al., 2004 a\ , the properties as a function of q of 
these regions coincide with those for the FA models, leading to the same estimates for 
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the scaling of the relaxation time. 

Relaxation processes in the Spiral model proceed in a broadly analogous manner. 
But the size of the critical defects that can expand further diverges already at q c , 
proportionally to the parallel length £11 of DP clusters. The typical size of incipient 
blocked clusters grows much more quickly, as L c ~ exp[c/(q — 9c) M ] with fi = uu (1 — z), 
where v ~ 1.73 is the critical exponent of £|| and z ~ 0.64 the exponent relating the 
parallel and transverse lengths of DP. As for FA models, since these clusters can be 
unblocked only from the boundary, r should diverge at least as L c . 

Interestingly, in the other model we have considered that has an ergodicity breaking 
transition at nonzero q c , the FA model on a Bethe lattice, simulation measurements of 
persistence and correlation functions ( |Sellitto et ai, 2005| ) show that relaxation times 
grow only as power laws on approaching the transition, r oc (q — q c )~ 7 with 7 ~ 2.9. 
Qualitatively, thus, the transition for FA models on Bethe lattices has the character- 
istics of a mode-coupling theory (MCT) ( |Gotze and Sjogren, 1992 1 arrest transition. 

Finally we discuss the East model, where the origin of the super- Arrhenius time 
scale is possibly easiest to see. The basic relaxation mechanism can be understood 
in terms of the dynamics of domains of occupied sites separated by empty sites. In 
the limit of small q one can then argue ( |Evans and Sollich, 1999[ ) that the typical 
relaxation time should scale as the minimal time required for an empty site to facilitate 
the motion of the first empty site to its right, which is typically at distance 1/q. The 
optimal path to create an empty site at a distance d involves an energy barrier of order 
log 2 (d) (the minimum over all paths of the maximal number of empty sites which we 
encounter along the path). Setting d= 1/q thus leads to the relaxation time estimate 
[/3 2 /ln2] ( |Evans and Sollich, 1999[ ). The above argument does 
not take into account the behavior on scales smaller than the typical distance between 
empty sites. Incorporating this turns out to halve the coefficient in the exponent, to 
t - exp [/3 2 / (2 In 2)] ( |Cancrini et ai, 2007| [Cancrini et al, 2008[ ). 

We have discussed so far only the overall timescale for the decay of correlation 
or persistence functions in KCMs, but not the time dependence of this decay. For 
the cooperative models one generically finds stretched exponential forms, as seen ex- 
perimentally, while non-cooperative models can exhibit power law tails reflecting the 
diffusive motion of the mobile defects. We refer to ( |Ritort a"n d Sollic h7"2003[ ) for an 
overview and quote only the example of the East model, where one can show that the 
stretching becomes extremely strong at low temperatures, with correlations and per- 
sistence both decaying as scaling functions of (i/r) 1 /^ 1 " 2 ! ( Evans and Sollich, 1999 
Buhot and Garrahan, 2001] |Sollich and Evans, 2003[ ). 



1.4.2 Some rigorous results 

In the recent years KCMs have been also analysed in the mathematical community. We 
summarize here those rigorous results which help to further understand the slow relax- 
ation of KCMs, and have in some cases corrected conjectures resulting from numerical 
simulations or intuitive arguments. 

The analysis of the large time behavior in the ergodic regime was started by 
( |Aldous and Diaconis, 2002] ) where for the East model it was established that [1/ ln2 — 
o(l)]ln(l/q) < ln(l/gap(g)) < [1/(2 In 2) + o(l)] ln 2 (l/g) with l/gap(g) the inverse 
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of the spectral gap of the Liouvillian operator generating the dynamics. The latter, 
which in a finite system is just the inverse of the smallest nonzero eigenvalue of the 
transition matrix, represents the longest relaxation time for all one-time quantities 
and so can be identified with the relaxation time scale r discussed above. The bounds 
of flAldous and Diaconis, 2002] ) quoted above then say that lnr scales for small q as 
ln 2 (l/<7) = j3 2 , with a prefactor between l/(21n2) and l/ln2, the upper bound be- 
ing the naive estimate of ( |Evans and Sollich, 1999[ ). These bounds were sharpened 
in (Cancrini et al, 2008), establishing that it is in fact the lower bound that gives 
the correct asymptotics: written in terms of the gap, liniq_>o hi(l/gap)/ln (1/q) = 
1/(2 In 2). 

Positivity of the spectral gap guarantees in particular exponential convergence 
in the large time limit: for any g(n) one has (g(n(t))g(n(0))) — (g(n(0))) 2 < const x 
exp(— 2gapi), where (. . .) is the mean over the initial Boltzmann equilibrium at empty 
site density q and over the stochastic process governing the evolution in time. In 
( |Cancrini et al, 2008[ ) it was shown that positivity of the spectral gap also guarantees 
exponential convergence of the persistence function P(t), which is the probability that 
a site does not change its state during a time interval of length t: P(t) < exp(— q gap t)+ 
exp(— (1 — g)gapi). In the same paper, a multi-scale approach was developed which 
allows one to prove positivity of the spectral gap in the whole ergodic region q > q c for 
all the choices of constraints described in Section fOl With this technique one can also 
derive flCancrini et al, 2008[ ) the following (sometimes optimal) bounds when q I q c . 
For FA-1 in d = 1, r oc 1/q 3 ; in d = 2, 1/q 2 < r < \n(l/q)/q 2 ; and in d > 3, 
l/q 1+2 / d < r < 1/q 2 . These bounds arc in agreement with the analytical results in 
(Jack et al., 200661, and confirm that the initial findings (Berthier et al., 200561 in 
d = 2 and d = 3, based on a mapping to DP, were incorrect. For the cooperative 
KCMs FA-2 and FA-3, results in ( |Cancrini et al, 2008| show that exp(q _1 ) < r < 
cxp(q~ 2 ) and exp[cxp(gf -1 )] < r < exp[exp(<7 -2 )], respectively, thus establishing a 
super- Arrhenius scaling compatible with QS Butler, 1991| [Reiter, 199"T| ). 

For conservative KCMs with non-cooperative behaviour, the diffusive scaling 1 / gap = 
0(L 2 ) in a volume of linear size L and the positivity of the self-diffusion coefficient 
at any density was established in (Bertini and Toninelli, 2004). Moreover the hydro- 
dynamic limit was studied there for a special class of models leading to a porous 
medium equation, namely a degenerate partial differential equation dtp = V(-D(p)Vp) 
with a diffusion coefficient D(p) vanishing as a power law of 1 — p when p 1. 
The methods used to establish these results use as a key ingredient the existence 
of a path between configurations which allows any two particles to be exchanged, by 
exploiting the presence of mobile macrovancancies. This approach, then, cannot be ex- 
tended to cooperative models. However, more recently a different technique has been 
devised ( |Cancrini et al, 2010a[ ) that proves also for cooperative models the diffusive 
scaling 1/gap = 0{L 2 ) and establishes in d = 2 the diffusive decay in time ~ 1/t of the 
density-density autocorrelation function. The self-diffusion coefficient for a specific co- 
operative model, namely the KA-m, was analysed in ( jToninelli and Biroli, 2004] ) where 
its positivity at any q was proved, modulo a conjecture on the behavior of a random 
walk in a random environment. 
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1.4.3 Non-equilibrium behaviour 

Due to space constraints we focus in this paper almost exclusively on the equilib- 
rium dynamics of KCMs. Non-equilibrium behaviour results if, for example, a non- 
conservative KCM is prepared in a configuration with a high density q of empty 
sites, and then q is reduced quickly. This corresponds to a sudden lowering of the 
temperature and so mimicks quench experiments in real glasses. If the final q 
is low enough, i.e. sufficiently close to q c , the time scale r for relaxation to the 
new equilibrium will be very long, and aging will occur: the properties of the sys- 
tem depend on the time elapsed since the quench. Aging can be monitored via two- 
time response and correlation functions, and one ask for example whether these are 
linked by a non-equilibrium fluctuation-dissipation theorem with some effective tem- 
perature ( |Crisanti and Ritort, 2003 ). Surprisingly, in some simple KCMs such as the 



East and FA-1 models, aging at low q is in fact easier to analyse than the dynamics 
at equilibrium flEvans and Sollich, 1 999 Ma yer and Sollich, 2007[ ). Work in this area 



up to late 2001 is summarized in (Ritort and Sollich, 2003), with more recent stud 



ies revealing a very rich phenomenology in the aging of KCMs ( |Mayer et al, 2006 



Leonard et al, 2007 Corberi and Cugliandolo, 2009). There has also been progress in 



rigorous approaches, with e.g. Refs. ( Faggionato et q/.niCancrmi et al, 2010&D estab- 



lishing the non-equilibrium behavior of the East model derived in ( Evans and Sollich, 1 999) 



1.5 Dynamical heterogeneity and its consequences 

The main success of KCMs has been the ability to account, at least qualitatively 
and sometimes also quantitatively, for many aspects of dynamical heterogeneity (DH) 
( |Ediger, 20 00; Glotzer, 2000] [Andersen, 2005[ ) observed in glass forming systems. DH 
emerges naturally in KCMs at high density or low temperature (Garrahan and Chandler, 2003; 
Pan et al, 2005[ ), and this is the topic of this section. We focus on the simpler KCMs, 



such as the FA-1, East and TLG constrained lattice gases. For these models results are 
easier to obtain both analytically and computationally than for, say, the Spiral model 
or the FA-m model with m > 1; nevertheless they display all the important physics. 
In particular, the dynamical scaling properties of the East model are close to those 



observed in realistic glass formers (see e.g. (Chandler and Garrahan, 2010 )) . 

The simplicity of KCMs allows one to make detailed predictions about the growth 
of dynamical correlation lengths and their scaling relations to growing relaxation times 



dGarrahan and Chandler, 2002||Toninelli et aZ.~ 004&; Whitcl anTef al, 2004|[Pan et al, 2005} 



|Toninelli et al, 200 5 b: Berthier and Garrahan, 2005 Jack et al, 2006 5| . These dynam 



ical lengthscales are an indication of spatial correlations that build up over time, 
so they are most easily extracted from multi-point correlation functions, which we 
discuss in Subsection 11.5.11 DH also tells us that the dynamics of glass formers is 
fluctuation dominated. A central consequence is transport decoupling, the break- 
down in standard transport relations of liquid state theory that are obtained un- 
der the assumption of homogeneous dynamics, such as that of Stokcs-Einstcin re- 



lating self-diffusion rate to viscosity (Swallen et al, 2003). KCMs provide a direct 



explanation for decoupling based on the response of a moving molecule to the distri- 
bution of local relaxation timescales in the host fluid flJung et al, 2004] ). Such mo- 
tion can be approximately quantified within a continuous-time random walk (CTRW) 
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formalism QBerthier and Garrahan, 2005[ ). Decoupling and CTRW are discussed in Sub- 
section [T32I Finally, the observed DH is only a mesoscopic phenomenon: dynamical 
lcngthscales are always finite (and transient) at non-zero temperature or at less than 
maximal density, that is, for q > q c (note that all models discussed in this section have 
q c = 0). Nevertheless, it is a precursor to a fully- fledged non-equilibrium, or "space- 
time", phase transition ( |Merolle et al., 2005||Jack et at, 200 6a; Gar rahan et ai, 2007[ ), 
which wc discuss briefly in Subsection 11.5.31 

1.5.1 Multipoint-correlations and susceptibilities 




Fig. 1.3 Dynamical Heterogeneity in KCMs. (a) Spatial distribution of the local persistence 
in the NEF model, at time ti/2 such that P(t\/2) = 1/2 (i.e., 50 of sites, shown in black, 
have flipped by time ti/2) at temperatures T — 1.0 (top) and T — 0.15 (bottom) for a 
system of size N — 40 3 . Adapted from Ref. QBerthier and G arrahan, 2005]), with permission, 
(b) Same, but for (2)-TLG model, at fixed density p — 0.77, and varying observation times, 
t = 10 3 , 10 4 , 10 s , 10 6 . From Ref. ( |Pan et al., 2005} , with permission. 

Figure 11.31 illustrates DH in KCMs. A convenient observable to quantify local 
relaxation is the persistence field, Pi(t) = 0,1, where indicates that site i has 
changed its state at least once up to time t, and 1 otherwise. The ensemble average, 
P(t) = (pi(t)) = limjv-voo N^ 1 J2iPi(t)> i s the persistence function discussed above. 
Fig. ll.3f a) shows the persistence field Pi (£1/2) in a three-dimensional version of the East 
model, or NEF (for North-East- Front model) ( |Berthier and Garrahan, 2005 ), where 



£1/2 is the time at which P{t\/2) = 1/2, i.e., half the system has relaxed and half has 
not. We see from Fig. I1.3f a) that relaxation is heterogeneous: there is a clear spatial 
segregation of sites which have relaxed, pi = (colored black), from those that have 
not, pi = (colored white). Relaxation dynamics is spatially correlated. Fig. I1.3f a) 
shows the persistence field at two different values of q (two different temperatures). The 
spatial extension of dynamic correlations increases with decreasing temperature, and 
therefore with increasing relaxation time. Note that these dynamical correlations are 
unrelated to thermodynamic correlations since the equilibrium measure of the NEF is 
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trivial at all temperatures. Fig. ll.3f a) shows similar persistence field plots for a density 
conserving model, the two- vacancy assisted (2)-TLG, but at a fixed high density p and 
for different observation times t. This figure shows that DH is a transient effect, with 
dynamical correlations becoming maximal at some intermediate time t* (see below). 

The extent of the dynamical correlations evident from Figure 1X731 can be quantified 
by means of multipoint functions ( Garrahan and C handle r" 2002| |Lacevic et ai, 2003| 



|Toninelli et ai, 2005 5| |Chandler et ai, 2006 1 . Consider for example the "four-point" 
structure factor, 

S 4 (k,t) = AT-W(t)X; [(pi(t) Pj (t)) - P 2 W] e&<?'-V (1.2) 

i,3 

This is the Fourier transform of the spatial correlation function of the persistence field. 
The factor Af(t) is a convenient normalization factor. We adopt the choice Af(t) = 
\P(t) — P 2 {t)] , which make £4 equal to unity if all the nonzero contributions come 
from "self" terms i = j. Fig. 11.4( a) shows S4 for the NEF. This is the structure 
factor of the DH pictures of Fig. 11.31 It has the characteristic shape of that of a 
system correlated over finite distances. The zero wave vector limit of S4 gives the 
"four-point" susceptibility, Xi{t) = S&(k — > 0,t), which provides an estimate of the 
correlation volume of DH. Fig. I1.4f b) shows Xi{t) f° r the NEF model as a function 
of observation time t, at different temperatures. Several things are apparent. \i is 
non-monotonic in time, indicating that DH is a transient phenomenon. It peaks at 
around the relaxation time of the persistence function. The peak value increases with 
decreasing temperature: dynamical correlations increase with increasing relaxation 
time. 

Multipoint functions reveal the scaling properties of DH. The dynamic suscepti- 
bilities defined above have their peak, x% = Xi(t*)i a t times t* close to the relaxation 
times r of the corresponding persistence functions. Scaling is controlled by the distance 
to the dynamical critical point at zero concentration of excitations or vacancies, so we 
expect xt ~ 9 -7 dToninelli et ai, 2005 fe| |Chandler et al., 2006[ ). Furthermore, we ex- 
pect the four-point structure factor to behave as S4 ~ xtf(k£), where / is a scaling 
function, and £ ~ q~ v a correlation length for DH at times t = t* . These exponents - 
and indeed scaling forms - may vary from model to model. For the simpler KCMs they 
can be calculated analytically. Figure \TA\ a) shows the numerical collapse of S4 in the 
NEF for 7 — 1 and v = 1/3 (Bcrth ier and Garrahan, 2005[ ). For other models, such as 
the FA-1 these exponents can be calculated analytically ( |Toninelli et al, 200 55). For 
example, in d = 3 one finds r ~ £ 4 . This is most easily understood from the fact that 
the upper critical dimension of the model is d c = 2, so exponents are d- independent 
above this flJack et al, 20066| . But in d = 2, the characteristic lengthscale is the dis- 
tance between vacancies of density q, £ ~ q^ 1 ^ 2 , and the timescale is related to this 
by the vacancy diffusion constant q/2 (see Section H.4. 1[) . giving r ~ £ 2 /(g/2) ~ £ 4 . 

Four-point and similar dynamical susceptibilities measure the dynamical fluctua- 
tions of global observables. These observables probe relaxation over a certain length- 
scale, so it is important to note that the scaling properties of the corresponding suscep- 
tibilities depend on such lengthscalcs flChandler et ai, 2006[ ). In the examples above 
the observable was the persistence function, which probes relaxation on a lengthscale of 
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Fig. 1.4 Dynamical correlations and scaling, (a) DH structure factor, Si{k,t±/2) in the NEF 
model at various temperatures. The data collapses under the scaling S4 — > qSi and k — > k£ 
with £ ~ q^ 1 ^ 3 , suggesting the values 7 = 1 and v = 1/3 for the scaling exponents (see 
text), (b) Time dependence of four point susceptibility Xi(^) m the NEF. Dynamical fluctu- 
ations probed by this function are maximal at time t* . The peak susceptibility \1 = X^(t*) 
grows with decreasing temperature. The peak time is approximately the relaxation time of 
the system, t* ~ r, as extracted from the persistence function, (c) Scaling of DH correla- 
tion lengths to relaxation times: FA-1 models follow a simple scaling law ( |Jack et al., 2 0066) 
(r ~ £ 4 for the three-dimensional FA-1 shown in the figure), while for East-like models 
r ~ £ z , where the dynamical exponent z(q) increases with decreasing q (i.e. decreasing 
temperature) ( |Sollich and Evans, 1999) . Adapted from Ref. ( |Berthier and Garrahan,~2 005 ), 
with permission, (d) The behaviour of four-point susceptibilities depend on the wavelength 
k of the observable which is used to define them. ^4(^1^) is plotted with its self-part, i.e. 
the terms with i = j in its definition, subtracted. From Ref. ( |Chandler et al., 2006[ ), with 
permission. 

one lattice spacing. Let us now consider instead the following susceptibility for the case 
of conservative (conserved density) KCMs, Xi(k, t) = Iinifc'->o ^2ij($Fi(k, t)SFj(—k, t)) 
where i,j label the N particles in the system, the position of the ith particle at time t 
is n(t), and SFi(k,t) = e ^-[r t (t)-v t (o)] _ ^k-[r t (t)- r m}) The wave vector k' appear- 
ing here is the analogue of k in S^(k,t) above. A normalization factor Af(t) could be 
included as before but this turns out to make no qualitative difference. This \i mea- 
sures the system to system fluctuations of the self-intermediate scattering function of 
wave vector F(k,t). It thus probes structural relaxation at lengthscales comparable 
to 2-7r/fc. Figure lOT d) shows how Xi(k,t) changes in behavior as we go from large 
to small k in the (2)-TLG: as we probe larger lengthscales (smaller k) Xi peaks at 
later times, and the initial power law growth changes exponent (Chandle r et al, 2006 ). 
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Fig. 1.5 Persistence/exchange decoupling. (a,b) Trajectory of a probe particle diffusing in a 
KCM. We show the case of the one-dimensional FA-1 for ease of illustration; similar behaviour 
is observed in other KCMs. The probe particle can make a diffusive step only if allowed by the 
local excitations in the KCM. (Sites with m = are shown as grey in the background.) At high 
temperatures/low densities (T = 3) diffusion becomes Fickian after short times and distances. 
Plot (b) is on the same scale as (a), showing that at lower temperatures (T = 0.8) diffusion 
is intermittent and non-Fickian over much longer timescales and lengthscales. (c) Timeline 
of displacement events. The waiting time between events is termed local exchange time; it 
is the time measured to the next event with the knowledge of when the previous one took 
place. The waiting time until the next event from an arbitrarily chosen starting observation 
time is termed persistence time. When the time series of events in non-Poissonian, as a 
consequence of dynamical correlation in the KCMs, typical persistence and exchange times 
are different: exchange times are dominated by the clustering of events, while persistence times 
are determined by the long quiescent periods. (d,e) The decoupling between the distributions 
of exchange (earlier curves) and persistence (later curves) times becomes more pronounced 
the lower the temperature. Results are shown for the East model, with T = f in (d) and 
T = 0.5 in (e). Adapted from Refs. ( |Jung et al, 20 04 Jun g et al, 2005] ) , with permission. 

This change in behavior is related to the non-Fickian to Fickian crossover of particle 
diffusion (Bcr thier et al., 2005 a[ ), which we discuss in the next subsection. 

1.5.2 Transport decoupling 

A central consequence of DH is transport decoupling. A prominent example is Stokes- 
Einstein breakdown flSwallen et al., 2003| : in deeply supercooled liquids the rate for 
self-diffusion is orders of magnitude larger than what would be predicted from the 
Stokes-Einstein relation between self-diffusion constant and viscosity, D s oc Sim- 
ilar transport relations of liquid state theory also break down near the glass transition 



(Chang and Sillescu, 1997 Ediger, 2000). This is a consequence of the dynamical fluc- 
tuations associated with DH. 

A major success of KCMs is the ability to rationalize this phenomenon. It does 
so in terms of the decoupling between the different fundamental timescales for local 
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relaxation ( [Jung et al, 2004[ |Jung et al., 2005[ ). How this comes about is illustrated 
in Fig. 11.51 Panels (a) and (b) show trajectories of a probe particle embedded in a 
KCM. This can be thought of as a molecule in a liquid that has been labelled for 
tracking while coarse-graininig over the rest of the system, thus describing its effective 
dynamics via a KCM. The motion of the probe particle is determined by the underlying 
fluctuations of the host KCM to which it is coupled. A natural dynamical rule is that 
the probe can make a diffusive jump from site i to site j only if both sites i and j are 
excited, n, = rij = ( [Jung et al, 2004[ ). Figs. 11.5( a) and (b) show the difference in the 
probe motion between high and low temperatures. At high T, Fig. I1.5f a). excitations 
in the KCM are plentiful and probe motion appears Brownian. At low T, Fig. II. 5f b), 
excitations are scarce, the dynamics of the KCM is heterogeneous, and probe motion 
is intermittent: the probe is immobile if immersed in an inactive space-time "bubble" ; 
in order to move it has to wait for an excitation to come along. 

There are two fundamental timescales that control this intermittent motion ( |Jung et al., 2004} 
Jung et al., 2005 1. The first one is the "persistence time" , t p , that is, the time the probe 



needs to wait to start moving for the first time, given an arbitrary start time for ob- 
servation. The second timcscalc is the (local) "exchange time", t x , the time between 
moves. Due to DH jump events are not Poissonian, but display "bunching" , and typ- 
ical persistence times can become much larger than typical exchange times, see Fig. 
1.51 c). This decoupling between persistence and exchange becomes more pronounced 
as temperatures is decreased, as illustrated in Figs. I1.5f d) and (e). 

We can approximately quantify the motion of the probe particle (Bcrthicr et al., 2005 a\ 



by means of a continuous-time random walk (CTRW) approach (Montro ll and Weiss, 1 965) 
The probe makes random walk steps of unit size at random times determined by the 
fluctuations of the host KCM. Lets assume that this random clock ticks at times drawn 
from the exchange time distribution, <j>(t x ), but which other than that are indepen- 
dently distributed. The probability for the probe to be at position r at time t, or 
van Hove function, is G* s (r, t) = X^m— o ^ m 

(t)T^(r) (we assume al = 1 for simplicity, 
extension to higher dimensions is straightforward). Here 7r m (£) is the probability that 
the probe made m steps after time t, and r^ m '(r) is the probability that a random 
walker is at a distance r after m steps. After Laplace transforming in time and Fourier 
transforming in space we obtain, 

F s (k, a) = P(a) + cos (k)M ' - fe> , (1.3) 

c 1 - cos (fc)0(er) 

where F s (k,o~) is the Laplace transform of the self-intermediate scattering function 
F s (k,t), which in turn is the Fourier transform of G s (r,t). Equation (|1.3[) is the 
Montroll- Weiss equation for the motion of the probe in the CTRW approximation. 
4>(o~), p(a) and P{o~) are the Laplace transforms of the exchange time distribution, the 
persistence time distribution, and the persistence function P(t), respectively. The last 
two functions appear in Eq. f| 1 . 3[) because the first step is determined by the persis- 
tence time: its distribution p(t p ) is related to the persistence function P(t) = no(t) by 
P(t) = J t °° p(t p )dt p . If the dyn amies is stationary, then there is nothing special about 
time zero and p(t p ) = (t x ) _1 4>{t')dt' , corresponding to a uniform average over all 
earlier jump times; the normalization factor is the inverse of the average exchange 
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time (t x ) = J °° dt(j)(t)t. From (JOJ) we can define the wavelength dependent relaxation 
time: r(k,T) = lim a ^o F s (k, a). We obtain ( |Berthier et al., 2005 a[ ) 

cos (fc) 



T(k) 



1 — cos (k) 



(1.4) 



where r x and r p are the average exchange and persistence times, respectively. At low 
temperatures we have persistence/exchange decoupling, t p 3> r x ( [Jung et al, 2005| ). 
The structural or alpha relaxation time is often defined as r Q = r(fe = 7r/a), where a 
is the lattice spacing which we set to unity. For these large wave vectors the first term 
of (| 1.4)1 dominates, and the structural relaxation time is set by the persistence time, 
T a ~ T p- For small enough wave vectors the second term in (|1.4I) dominates, r(fc) ss 
r x /k 2 , and since the limit of k —¥ defines the diffusion rate we find that D sa t" 1 . 
This explains Stokcs-Einstein breakdown at low temperatures: Dr a w t p /t x ^ const 




(b) 



10 10" 10 10 10 




Fig. 1.6 Transport decoupling in KCMs. (a) Scaling of probe or self-diffusion constant with 
structural relaxation time in the one-dimensional FA-1 and East models (left scale) and 
two- and three- dimensional lattice gases (right scale). Solid lines indicate power law fits 



low enough q, indicating a breakdown of the Stokes-Einstein relation, as a consequence of 
dynamical fluctuations. (For the FA-1 model 5 — 1 for dimensions larger than its critical di- 
mension d c = 2 ( |Jack et al, 20066} .) Adapted from Refs. < )Jung et al., 2004|[Pan"ef al, 2005) 
|Ashton, 2009} , with permission, (b) Lengthscale dependent relaxation time r(fc) in the FA-1 
and East models. At short lengthscales the relaxation time is k independent and determined 
by the persistence time, r(fc) ~ r p . At larger lengthscales it becomes diffusive, r(k) ~ r x /fc 2 . 
This crossover is controlled by the Fickian lengthscale I* ~ ^/r p /r x . The figure shows that 
r(k) at different q collapse under k — ¥ kl* . From Ref. ( [Berthier et al., 2005 a\ , with permission. 

The decoupling between persistence and exchange times is an effect of dynamical 
fluctuations. At low temperatures the self-diffusion constant seems to scale as a frac- 
tional power of the relaxation time D ~ t~ s , with S < 1 (the Stokes-Einstein relation 
is 5 = 1). This is the case for all KCMs with strong enough constraints, as shown 
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in Fig. 11.6( a) for the FA-1 in d = 1, the East model and various constrained lat- 
tice gases ( |Jung et ai, 2004[ |Pan et ai, 2005| |Ashton, 2009[ ) . The fractional exponent 
S ~ 0.6 — 0.8 is not distinct from that observed in experiments ( |Swallen et ai, 2003[ 



Swallen et al, 2009 ). 

Equation (|1.4[) indicates that there is a spectrum of timescales that interpolate 
between a k independent value r p at shorter lengthscalcs, to a diffusive timcscalc 
T x /k 2 at large lengthscalcs ( |Berthier et al, 2005 a| ). This crossover is shown in Fig. 
1.6( b) for the FA-1 and East models. It is the crossover from non-Fickian diffusion at 



short lcngthscales to eventual Fickian diffusion at long enough ones. The length I* at 
which this crossover takes place, or "Fickian lengthscale" , is given by /* cx \/t p /t k , 
and grows with decreasing temperature/increasing density. It is the distance a particle 
has to move before it forgets how long it took to make the first step. Fig. 11.6( b) 
shows how r(fc) at different temperatures collapse under k — > l*k. The non-Fickian 
to Fickian crossover is also responsible for the wavelength dependence of four-point 



functions flChandler et al, 2006[ ), Fig.[Tl|d). 

The CTRW analysis can be extended to describe the effect of driving, for example 
by externally forcing the probes. The competition between timescales in this case 
leads to interesting non-linear response behaviour, such as non-monotonic differential 
mobility and giant diffusivity QJack et al, 2008[ ). Furthermore, a study of (|1.3p in the 



crossover regime between non-Fickian and Fickian explains (Chaudhuri et al, 2007) 



the exponential tails observed (Stariolo and Fabricius, 2006) in van Hove functions at 



intermediate times. The waiting time distributions used in the CTRW analysis above 
are the ones that are obtained from the study of KCMs. The CTRW approach can 
also be used by assuming a different origin for the waiting time distributions, as for 
example in the analysis of metabasin transitions ( |Heuer, 2008[ ) in atomistic liquids. 

1.5.3 Space-time phase transitions 

Thcrmodynamically KCMs are trivial, so all interesting behaviour is dynamical. Nev- 
ertheless, DH pictures such as those of Fig. 11.31 are suggestive of phase separation 
between two distinct phases. The phases of Fig. 11.31 are distinguished by their dy- 
namics: dark regions are dynamically active while light ones are dynamically inactive. 
Furthermore, the phase separation in the spatial projection of low temperature/high 
density equilibrium trajectories, such as those of Fig. 11.31 is only mesoscopic: when 
coarse-grained over large enough lengthscales, irrespective of time, the space projected 
trajectories are homogeneous. Here we show how these observations are directly related 



to a true non-equilibrium phase transition (Garrahan et al, 2007), which in contrast 
to thermodynamic transitions, occurs in ensembles of trajectories and is driven by 
non-equilibrium driving fields. This phase transitions can be studied by recourse to 
the large-deviation method ( |Lecomte et ai, 2007[ |Touchette, 2009[ ). 

A convenient order parameter to discern active and inactive dynamics is the "dy- 
namical activity" K ( jGarrahan et ai, 2007P , defined as the total number of configu- 
ration changes in a trajectory. In a non-conservative KCM it would amount to the 
total number of local changes from empty to occupied or vice versa, and in a lattice 
gas (conservative KCM) to the total number of particle displacements. The activity 
is extensive in space-time volume, i.e., typically K = 0(Nt), where N is the number 
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(a) (b) (c) 




Fig. 1.7 Space-time phase transitions in KCMs. (a) Variational dynamical free-energy J-(n) 
in terms of mean excitation density. Kinetic constraints give rise to a bistable Tin) (full 
line). In the absence of kinetic constraints the corresponding variational function is unistable 
(dashed line), (b) Mean-field estimate of large-deviation function tp{s). There is a singularity 
at s = 0, indicating a dynamical phase-transition, (c) Mean activity K(s) = —ij)'(s) as a 
function of s. This dynamical order parameter shows a discontinuous jump at s = 0: the 
transition is a first-order one between an active dynamical phase and an inactive dynamical 
phase, (d) The dynamical first-order scenario is also present in finite dimensions, as shown 
from numerical computation of ip(s) and K(s) for various KCMs. (e) Dynamical phase di- 
agram. The line s = is one of first-order coexistence between the active dynamical phase 
(s < 0) and the inactive one (s > 0). It extends all the way along the q axis. The critical 
point at q = controls the scaling behaviour discussed in Subsection 11.41 Adapted from 
Refs. (fGarraha n et al., 2007[ [Garrahan et al, 2009 a\ , with permission. 

of lattice sites and t the time extension of the trajectory. Each trajectory x(t) in the 
ensemble of (equilibrium) trajectories of length t has a total activity, K[x(t)]. The ac- 
tivity is thus distributed, Pt(K) = (S (K - K[x(t)]\\ where the average is over the 
set of equilibrium trajectories, {x(t)}. At long enough times this probability acquires a 
large-deviation form flLecomte et al, 2007| |Touchette, 2009] ). P t (K) w e -«v(*7*). The 
function ip(k) is called a large-deviation function, and it plays in this dynamical context 
the same role as, for example, the entropy density in the micro canonical ensemble of 
equilibrium statistical mechanics. Alternatively we can consider the generating func- 
tion of K, 

K 

which also displays a large-deviation form. The large-deviation function tp(s) is akin 
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to a free-energy density, and is related to ip(K/t) by a Legendre transform, ip( s ) = 
— minfe [ip(k) + sk]. Just like a free-energy in a thermodynamic problem, the function 
ip(s) carries the information of dynamical phase behaviour. Specifically, its singularities 
indicate dynamical phase transitions ( |Garrahan et al., 2007[ ) . 

The calculation of the large-deviation function ij)(s) is simplified greatly by the 
following observation QLebowitz and Spohn, 1999] ) : if W is the master operator that 
generates the stochastic dynamics, then ip(s) is the largest eigenvalue of a modified 
operator W s , where Wo = W. This reduces the calculation from that of computing a 
"partition sum" , Eq. (|1.5[) , to an eigenvalue problem. For example, for a spin facilitated 
model, W s = £\ fi(n) {e~ s [(1 - p)af + po~] - pn { - (1 - p)}, where fi(n) is the 
kinetic constraint on site i, and of are the raising/lowering operators on site i. While it 
is not always possible to diagonalise such an operator analytically, bounds for its largest 
eigenvalue can be estimated variationally. For the FA-m this amounts to minimising a 
Landau free energy F(n) = —n m (2e~ S y/qn — q — n). The factor n m comes from the 
kinetic constraint, and makes F(n) non-linear enough to allow for multiple minima, 
Fig. 11.7( a). The corresponding ip(s) has a singular structure, the first derivative being 
discontinuous at s = 0, Figs. 11.7( b) and (c). The meaning of this is the following. 
Trajectories are organized into two dynamical phases, an active one with K > and an 
inactive one with K — 0. The field s determines the bias for or against activity. For s < 
the active phase is the dominant one, while for s > the inactive phase dominates. At 
s = the probability of trajectories in either phase is equal in the t — > oo limit and we 
have dynamical first-order phase coexistence ( |Garrahan et al., 2007[ ). This situation 
occurs in all KCMs, Fig. 11.7( d). Interestingly, similar dynamical phase structure is 
observed in spin glass mo dels (|van Duijvendijk et al., 2010 Jack and Garrahan, 2010 ) 



and in atomistic liquids (Hedges et al., 2009 ) 



Actual dynamics takes place at s = 0. The results above show that this is the 
condition for dynamical coexistence in the bulk, i.e., infinitely far from boundaries 
active (crgodic) or inactive (non-ergodic) trajectories are equally likely. However, just 
like in the case of ordinary phase transitions, boundary fields can bias the bulk into 
one of the coexisting phases. In the case of dynamics, initial conditions play the role 
of a (time) boundary. In particular, almost all possible initial configurations chosen 
from the equilibrium static distribution at non-zero temperature or less than maximal 
density will select the active dynamical phase. In this case the inactive phase manifests 
via rare region effects, giving rise to DH only at the mesoscopic scale. An important 
open question is whether there are physical controllable fields that play the role of s 
in the analysis above. 

1.6 Summary and outlook 

In this chapter we have attempted to summarise recent developments in the study of 
KCMs as models of glass formers. In the long tradition of statistical mechanics KCMs 
provide simplified models that capture important ideas about the fundamental physics 
behind the phenomenology of glassy systems. Their simplicity allows for detailed study, 
which in turn gives rise to further physical insights into the glass transition problem. 
The central message from KCMs is that the complex and cooperative dynamics of 
glass forming systems can be achieved without recourse to complex thermodynamic 
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behaviour: in KCMs thermodynamics plays essentially no role, and complex dynamics 
emerge from rather simple local kinetic rules. These rules are local and free of disorder, 
but nevertheless give rise to dynamical frustration. 

The irrelevance of thermodynamics for glassy dynamics that the study of KCMs 
suggest contrasts sharply with approaches such as that of the random first-other tran- 
sition theory ( |Lubchenko and Wolynes, 2007| |Mezard and Parisi, 2000[ ) where ther- 
modynamics is essential. Whether thermodynamic aspects are relevant or not to glass 
transition phenomena is still a matter of debate, but to the extent that they are KCMs 
can say very little about them. This can either be seen as a flaw of the KCM based 
approach ( |Biroli et al., 2005"] ) or as an indication that these aspects are described by 
degrees of freedom that do not contribute too much to the long time dynamics and 
have therefore been coarse-grained out ( |Chandler and Garrahan, 2005[ ). Furthermore, 
by their coarse-grained and lattice based nature KCMs can in principle say very lit- 
tle about short distance/short time dynamics, such as beta- relaxation or anomalous 
vibrations. There is however evidence that these short scale phenomena are coupled 
to longer scale dynamic heterogeneity (see for example ( |Widmer-Cooper et al., 200"H 
|Brito and Wyart, 2007[ )) so it may be possible to capture some of these effects with 
generalisations of KCMs ( |Moreno and Colmenero, 2007| |Ashton and Garrahan, 2009[ ) . 
In any case, KCMs provide an explicitly real-space picture of glassy dynamics. Their 
main success has been the rationalization of dynamic heterogeneity. While DH can 
be analysed with other approaches, such as generalisations of mode-coupling theory 
( |Biroli et al., 2006] ), the immediacy of the results and explanations for DH related 
phenomena obtained from KCMs is remarkable. 

While we know a lot about KCMs we still do not have a satisfactory under- 
standing of how they emerge as an effective description from realistic systems. (This 
is also the case in alternative approaches, be it the random first-order transition 
( |Lubchenko and Wolynes, 2007| , or frustrated limited domains (Kivclso n and Tarjus, 2008] ), 
where the idealised models that display the proposed behaviour cannot be readily ob- 
tained from realistic liquid systems.) It is usually argued ( Garrahan and Chandler, 2003 1 
that KCMs ought to emerge from some form of local coarse-graining of a microscopic 
system, but this procedure has not been shown to work just yet (except in highly sim- 
plified situations (Garraha n and Newman, 2000| Garrah an et al., 2009 fc[ )). Proving a 
direct connection between atomistic liquids and KCMs is arguably the central open 
problem in this field. 
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